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I. INTRODUCTION 

In recent years there has been a continuous interest in the study of excited heavy 
mesons composed of a heavy quark and a hght anti-quark. This interest arises from 
several reasons. Some of these excited states have been observed in experiments. They 
will be objects of further study in future experiments with B-factories. In particular, they 
are useful for tagging in CP experiments [0]. Theoretically, the relative simplicity of the 
dynamical condition makes it of interest for exploring the internal dynamics of systems 
containing a light quark. 

Important progresses in the theoretical description of such systems have been achieved 
with the development of the Heavy Quark Effective Theory (HQET) 0. Based on the 
spin-flavor symmetry of QCD, valid in the limit of infinite heavy quark mass tuq, this 
framework provides a systematic expansion of heavy hadron spectra and transition am- 
plitudes in terms of the leading contribution, plus corrections decreasing as powers of 
l/rriQ. However, to obtain some detailed predictions one needs to combine it with some 
non-perturbative methods. The spectra and decay widths of heavy meson excited states 
have been studied with the l/rriQ expansion in the relativistic Bethe-Salpeter equations 
in P,^. They can also be studied with the QCD sum rules in HQET which have been 
used both for ground states of heavy mesons ^^ and for lowest excited heavy meson 
doublets (0+, 1+) and (l+,2+) at leading order of the l/mg expansion [§. For related 
results, see 0. 

The effective Lagrangian of the HQET, up to order l/rriQ, can be written as 

C,s = KivDK + -^IC + -^S + 0{l/ml), (1) 

ZrriQ Ztuq 

where /i^(x) is the velocity-dependent field related to the original heavy-quark field Q{x) 
by 

^.(x)=e^™«"-^^^Q(a;), (2) 

/C is the operator of nonrelativistic kinetic energy with a negative sign defined as 

}C = K {iD^fK , (3) 

where D'^ = D^ — {v ■ D) v^^ with D^ = d^ — ig A^ is the gauge- covariant derivative, and 
S is the Pauli term, describing the chromomagnetic interaction: 



5 = - Cmag{mQ/fl) h^ (Xf^^G^^hy , (4) 



r \ \ 3//3o 



where Cmag = r-r- , /^o = H ~ 2nf/3. Apart from leading contribution, the 



Lagrangian density contains to OII/itlq) accuracy two additional operators /C and S. 
The matrix elements of these two operators over heavy meson states are fundamental 
parameters in the l/rriQ expansion with HQET that should be either determined phe- 
nomenologically from experimental data or estimated using a nonperturbative theoretical 
approach. The matrix element of the chromomagnetic interaction operator is the leading 
contribution to the mass splitting among the states belonging to the same heavy meson 
doublet. 

In our previous work §], the masses of lowest excited heavy meson doublets (0^, 1+) 
and (l+,2+) are calculated at the leading order of I/ttlq expansion by using the QCD 
sum rule approach in HQET. As an extension of this work, in the present article we shall 
use the QCD sum rules in HQET to obtain quantitative estimates of the matrix elements 
of the kinetic energy and chromomagnetic interaction operators, which are responsible for 
the kinetic energy and the chromomagnetic mass corrections. 

The remainder of this paper is organized as follows. In Section |l| we begin with a 
brief review on the interpolating currents for excited heavy mesons and some important 
properties of these currents. Some results of two-point functions for the doublets (0"^, 1"^) 
and (1^, 2"'') at the leading order of the l/rriQ expansion are outlined. Section |T| is 
devoted to 0{l/mQ) corrections, the sum rules for the relevant matrix elements of kinetic 
energy and chromomagnetic interaction operator are derived. Finally, numerical results 
are presented in Section |V[ 

II. INTERPOLATING CURRENTS FOR HEAVY MESONS OF ARBITRARY 
SPIN AND PARITY AND TWO-POINT CORRELATION FUNCTION 

A basic element in the application of QCD sum rules to excited heavy mesons is to 
choose a set of appropriate interpolating currents in terms of quark fields which create 
(annihilate) definite states of excited heavy mesons. For excited heavy mesons with 
arbitrary spin and parity this problem has been studied in our previous paper §]. Here 
we briefly outline the results of this analysis which will be used in this article. 



From the general form of the leading order B-S wave function obtained in [0, we can 
write the general expression of the interpolating current creating a excited heavy meson 
with arbitrary spin j and parity P as 



j'^ai-'-Oj 



4«i--aj}, 



j;:^r(^) = h.{x)Tj^:r'(T^.Mx) 



(5a) 



or 



j;-.-"^(a;) = Kix)T\'V,p\v^,)i-^) p,,q{x) 






(5b) 



where i = 1,2 correspond to two serieses of doublets of the spin-parity [j^ ^^^ , (j + 

d 

1)*^~^^^ ] and [j^~^^^ , {j+l)^~^^^] respectively, V^^ is the covariant derivative 7- igAt{x), 

OXt 

subscript t denotes the transverse component of a 4- vector and 

r{a,....,}(p^j = Symmetrize |r"i-"^(P,J - ^^r"'^l;4r"i"'2"3-.(2)^j| 

with gf^"^'^ = (y(°i"2 _yai^a2^ Pqj, ^Yie doublets of the spin-parity [j^~^^^ , {j + l)^~^^^ ] and 
[j^~^^\ (j + l)^""'^^''], the expressions for T°'^"'°'^ (T)^^) have been explicitly given in [3, P^] as 



r(^.J 



r(^.J 






' ^rp.j ioTj^-^y^' 



2j + l 



I2j + 1 



-tyv:^...v:^AK: 
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I V 2j + 2 ' ' ^' ^* ' ^* 2j + 1 

The currents hermitian conjugate to (^) and (|5bD 



ir V. 



for(j + l)(-i)^^' 

for (j + 1)(-^)' 
forj(-i)^ 



(6) 



(7) 



a::t> 



j;:^r^(x) = g(a:)fj;^V-"^^(5.J/..(.), 



(8a) 
(8b) 



where T = 7°r''"7°, correspond to the annihilation of the same meson. 

The currents defined in (|^) and (|^) have nice properties. We had shown in [^] that 
two currents in the sets (|5a|) and (pE|) with not identical values of j, P, i never mix 
in the mg -^ 00 limit either in the heavy meson or in the quark-gluon level and heavy 
quark symmetry is explicit with these currents. This verifies that they are the appropriate 



interpolating currents for heavy meson states with definite j, P and the hght quark angular 
momentum ji which is conserved in QCD in rriQ -^ oo limit. 

The properties for currents mentioned above are important for applications to QCD 
sum rules for excited heavy mesons. If one use other currents or study the sum rules 
in full QCD, there are in general contributions from two nearby poles corresponding to 
states of the same j, P which have different values of the total angular momentum of 
the light component ji = j ± - in the rriQ -^ oo limit. Their contributions may not be 
separated correctly. Furthermore, the mixing of such two states can be calculated within 
our formalism by introducing the Oll/niQ) terms in the Lagrangian in HQET. 

Let \j, P, i) be a heavy meson state with the quantum numbers j, P, i in the rriQ -^ oo 
limit, together with the corresponding interpolating current, we have [§] 

{0\j;yp{0)\j\P\t') = fp,,5,,>5pp>5u't'-"'^ , (9) 

where the decay constant fjpi = fpj^ has the same value for the two states in the same 
doublet and r]"^""^ is the transverse, symmetric and traceless polarization tensor. 

In the framework of QCD sum rules the decay constant can be obtained from the 
two-point correlator: 

^^■a„P,-P, (fc) ^ , I rf4^g.fc.x^Q|y ^^i^p {X)J]^-'^' (0)) |0) , (10) 

where k is the residual momentum, corresponding to the decomposition p^ = rnqv^ + k^ 
of the heavy meson momentum. This expression has the following form at the leading 
order [^ 

KpP''''"'m = {-P S gt^'^ ■ --gT'^ n,,p,(o;) , (11) 

where cu = 2f ■ /c is twice the external off-shell energy and S denotes symmetrization and 
subtracting the trace terms in the sets {a\ ■ ■ -aj) and (/^i ■ ■ ■ fij)-, and Tl^^p^iiuS) = Upj^^u) 
is independent of j and satisfies dispersion relations of the form: 

IljP^i(iu) = / du •'' '^ h substractions , (12) 

where pjpp) = — Imn,-pj(z/) represents the spectral density. 

n 

We shall confine our study to the doublets (O"*", l"*") and (1"'",2"'") here. According to 
(^), (|5bD and (|^), there are two possible choices for currents creating O"*" and 1~*" of the 
doublet (0+,l+), either 



or 



Tt 

-'0,+ 


2 — 


^M, 


Tt" 


2 ^ 


^ /i.7'7rg , 


/'t _ 


1 


K{-i) ptq , 


/t _ 


1 

V2 


/i,7'7r(-^) A? • 



(13) 
(14) 



(15) 
(16) 

Similarly, there are two possible choices for the currents creating 1+ and 2^^ of the doublet 
(l+,2+). One is Q 

4,+,i = sjl h.l'H) (Vt - ^7? Pt) q, (17) 

ft"!. "2 _ . / 7, V J_ I ^,ctiT^a2 _|_ ^,a2T^ai ^ „aia2 



4+:r = \hK^ ( 7rA"^ + 7^^"^ - 3 ar^ pt)q- 



Another choice is obtained by adding a factor —i pt to (p^Tf ) and (0). Note that, without 
the last term in the bracket in ( pTj) the current would couple also to the 1+ state in the 
doublet (O"*", l"*") even in the limit of infinite tuq. 

Usually, the currents with the least number of derivatives are used in the QCD sum 
rule approach. The sum rules with them have better convergence in the high energy region 
and often have better stability. However, there is a motivation for using the currents (pISj), 



(P^) for the doublet (0^,1^). In the non-relativistic quark model, which usually gives 
correct ordering of energy levels of hadron states, the doublets (0"^, 1^) and (1+, 2+) are 
orbital p-wave states which correspond to one derivative in the space wave functions. In 
fact, as proved in 0, the coupling constant for the currents (|15|), (|T6|) and that for (p!7D, 
(|l^) are proportional to the large components of the B-S wave functions for the doublet 
(O'^, 1^) and (1^, 2"'') respectively in the non-relativistic approximation. Therefore, we 
shall consider both the currents (|T3]), (^^ and (|TB|), (|TB|) for the doublet (0+, 1+). 

For the doublet (0^,1^*"), when the currents J'0,+,2, J'1,+,2 in (HI), (111) are used the 



sum rule (same for the two states) after the Borel transformation is found to be 



^There is a typo in the coefficient of the last term of Eq. (25) in ref. [8], the coefficient 2 should 
be 2/3. 



/^e-^/^ = ^ f u'e-/^du -^^ml (gg) . (19) 



The corresponding formula when the current Jo,+,2 and Ji,+,2 in (HI) and ([1^) are used 
instead of J'0,+,2 and J'1,+,2 is the following 

/^e-^/^ = ^ f uj\-l^du + \ m - ^ m^ m . (20) 

When the currents (0) and ([T8|) are used the sum rule for the (l+,2+) doublet is found 
to be 

Here m-Q (gg) = (qgo'^^G'^'^q). In the above derivations, we have confined us to terms of 
the lowest order in perturbation and operators of dimension less than six. Above results 
shall be used in the next Section. 



III. THE SUM RULES AT THE 0{1/Mq) ORDER 

From Lorentz covariance, the correlator ( p!OD still has the form (|Tl|) beyond the leading 
order of l/mq expansion if we confine fc to be a longitudinal vector. But now Iljpjj{uj) 
depends on j due to the chromomagnetic interaction operator S. Inserting the heavy 
meson eigen-state of the Hamiltonian up to the order 0{l/mQ), the pole term on the 
hadron side becomes 

U(u) (/ + ^/)' f" 2^^/' I 2/^/ r22) 

^^{^)poie 2{A + 6m) -u 2k -u {2k-ujf^ 2k-uj' ^ ' 

where Sm and 6f are of the order 0{l/mQ). 

To extract 6m in ( ^2]) we follow the approach of and consider the three-point 
correlation functions 

SoU2pP'''-'^{uJ,uj') = i'' j d'xd'ye''-^''-ym{j;:p-p{x) 0(0) j]%:'^{y)) |0) , (23) 

where O = K, 01 S. The sum 5^11 + ^^II is equal to the correlator ([T0|) in the order 
0{l/m,Q) when uj = u . The scalar function corresponding to (p3|) can be represented as 
the double dispersion integral 



fonKJ) = l/ /°''-;;'f"'-' . (24) 

TX^ J [S — UJ){S' — UO ) 



The pole parts for them are 






PK 



(2A-c<j)(2A-cj') 



2A-UJ 2A-cj' ' 
Gsioo') , Gsiuj) 



duf 



2A-iu 2A - w' 



(25) 
(26) 



where 



2dM'>^j,p,ji = {j,P,Ji\K9(^f„yG^''K\j,P,ji) , 

djM = djjj, dj^_i/2,ji = 2ji + 2, dj^+i/2jj = -2ji. 



(27) 
(28) 
(29) 



Let u = uj in 



and (|26D, and compare it with ( ^2]) one obtains 
1 



(5m 



4:171, 



-{K + dMCmag^) ■ 



(30) 



The simple pole term in ( pSD and (p6D comes from the region in which s{s') = 2A and s{s') 
is at the pole for a radical excited state or in the continuum. These terms are suppressed 
by making double Borel transformation for both u and u'. One obtains thus the sum 
rules for K and S as 



J2^g-2A/T^ 



cjc r^< 



JO 



JQ 



dujduj'e-^^+^'^/^^ps{uj,uj') 



(31) 
(32) 



where the spectral densities are obtained from straitforward calculations in HQET. 

Confining us to the leading order of perturbation and the operators with dimension 
D < 5 in OPE, the relevant Feynman diagrams, which contribute to ^y^H and SgU are 
shown in Fig. 1 and Fig. 2 respectively. The graphs missing there turn out to be 
vanishing in the fixed point gauge x^A^{x) = [jlT|. The spectra functions p/c and ps 



can be calculated from these diagrams. For the jf = - doublet, we find 



/2^g-2A/T 



277r2 Jo 



u 



e-''/^du;+^{asGG)T 



2% 



fSe- 



-2A/T 



487r 



{asGG) T^ 



(33a) 
(33b) 



when the currents (|T^) and ([T6|) are used and 



/2^g-2A/T 



2%2 Jq 



cu^e-^/'^du 



25 



TC 



-{asGG) T - - m^ (qq) 



f^e-^-^/^ = ±^{a,GG)T + ^ml{qq) 



(34a) 
(34b) 



when the currents (|T3|) and (|T^ with no extra derivative are used. In contrast to ( P3b|) , 
there is a nonvanishing mixing condensate term in (|34b|) . It arises from the diagram 



depicted in Fig. 2(b). For jf = - doublet, a straitforward calculation yields 



/2^g-2A/T 



2^2 



LU^e-'^/^du; 
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3 X 2% 



{asGG) T^ 



fEe-^^/^ = -^{asGG)T\ 



727r 



(35a) 
(35b) 



Combining (^3]), (0) and (|35|) with (plQ]), (pO|) and (^) we can obtain sum rules for K 
and S in the three cases. 

The above results are not the whole story. The spin-symmetry violating term S not 
only causes splitting of masses in the same doublet, but also causes mixing of states with 
the same j, P but different ji. This mixing is characterized by the matrix element 



This quantity can be extracted from the correlator 



Jd'xd'ye''-^-'''y{0\T (j,,p,,+ i(x) K ^-cr^.G^^hM ^1p,,-(i/) ) |0) 



(36) 



(37) 



the double pole term of which is 



(2A_ 



cu)(2A_ 



UJ' 



(3J 



-J+2 '^ +'1-2 

With the methods similar to that used above one can obtain the sum rule of m^^lj^ 
For two 1"*" states we find 



As an example, here we have used the currents ([161) and (|l^) to do the analysis. 
The masses of two 1^ states are obtained by diagonalizing the mass matrix 



(39) 



A+i 



4mo +'2 " +'2 



Cmagm-i i(l" 



4mQ 

1 



Am. 



Q 



^'2 Am, 



G-niagl^klij' 



A-f^-l-^l + 5 Gmag^^ 



\ 



(40) 



The correction to the masses due to the mixing is formally of the order of 0(l/mQ). How- 
ever, due to the smallness of the difference between the two diagonal elements there is 
some possibility that this correction is numerically comparable to the 0{l/mQ) correc- 
tions. Therefore we shall calculate it also. 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

We now turn to the numerical evaluation of these sum rules. We can eliminate the 
explicit dependence of the sum rules obtained in Section |ITT| on / and A by means of 
the two-point function sum rules outlined in Section |I[ Dividing the sum rules in (|33D, 
(0) and (^) by the sum rules in (p!9|), (pO|) and ([21| ) respectively, we obtain relevant 
expressions for the K and S as functions of the Borel parameter T and the continuum 
threshold Uc- For the QCD parameters entering the theoretical expressions, we use the 
following standard values 

(qq) = -(0.24 GeV)^ , 
{asGG) = 0.038 GeV^ , 

ml = 0.8 GeV^ . (41) 

The results for K as a. function of the Borel variable T, and for different values of the 
continuum threshold Uc, are shown in Fig. 3. In particular, in Fig. 3a and 3b we depict 
the results of two kinds of sum rules for doublet (0^, 1^), which correspond to using two 
kinds of currents ( [13|) ([l^) and (|l^) (|16D respectively. In Fig. 3c, we display the result of 



sum rule associated with doublet (1+, 2+). 

Following the QCD sum rule procedure, we should check the existence of stable regions 
for the results of K in the Borel variable T and the continuum threshold Uc. As one can 
see from Fig. 3, stability starts at values of Borel parameter T slightly larger than 0.7, 
and stretches practically to T -^ oo. The familiar criterion that both the high-order 
power corrections and the contribution of the continuum model should not be too large, 
generally speaking, less than 30%, rather strongly restricts the 'working' region. Thus, 
the stability at large T values is not useful, since in this region the sum rules is strongly 
contaminated by higher resonance states. This continuum model contamination problem 
is quite severe also in the case of sum rules analysis of ground state heavy meson to I/uiq. 
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|P,0]. It originates from the high power dependence of the spectral densities. According to 
1^, it is better to determine the working window by considering the stabihty of two-point 
sum rules, and then to use the same values in the evaluation of the three-point sum rules. 
The choice of so-called stability window does not necessarily coincide with the stability 
plateau for the three-point sum rules. In the case of two-point function sum rules one 
finds that the working windows are Uc = 2.5 — 2.9 GeV, T = 0.7 — 0.9 GeV for the doublet 
(0+, 1+) and cje = 2.7 - 3.2 GeV, T = 0.7 - 0.9 GeV for the doublet (1+, 2+) §. 

From Fig. 3, taking the working regions of the sum rules as that for the two-pont 
function for the three cases, we obtain the set of values for K as following: 

K = -1.85 ±0.30 GeV\ (42) 

for the doublet (0+, 1+) when the currents ( pTSf ) (|T^ without the derivative are used, 

K = -1.90 ± 0.40 GeV^ (43) 

for the same doublet when the currents ([T5|), ([T6|) with the derivative are used, and 

K = -2.00 ± 0.40 GeV^ (44) 

for the doublet (l+,2+). 

The errors reflect the variations with the Borel parameter T and the continuum thresh- 
old uJc within the working windows. As can be seen when comparing Fig. 3a and Fig. 3b 
with each other, K value in Fig. 3b is more sensitive to the value of uJc than that in Fig. 
3a owing to the higher dimension of spectral density for the using of current with extra 
derivative. 

In Fig. 4 we show the numerical results of the sum rules for E, corresponding to the 
three chosen values of the continuum threshold Uc- From these figures, we see that there is 
no stability plateau for the sum rules of S. This is related to the absence of perturbation 
terms for S in our approximation. However, due to the strategy mentioned above we can 
take the working window as that of two-point function sum rules like the case of K. The 
results for S for the doublet (0+, 1^) in the working regions of the sum rules are 

E = 0.023 ± 0.003 GeV^ , (45) 

when the currents (|l^) (jl^) without the derivative are used and 

10 



S = 0.014 ±0.003 GeV^ , (46) 

when the currents (|T3p, (|T^) with the derivative are used. For the doublet (l+,2+), the 
value of S is given by 

S = 0.020 ± 0.003 GeV^ . (47) 

Dividing the sum rule in ( pQ]) by the sum rules in ( pOD and ([21|), we also obtain the 
expression for mi 3(1"'"). In Fig. 5 we plot the result of this sum rules as a function of 
T, and for three values of Uc- Using the working windows as cUc = 2.6 — 3.0 GeV and 
T = 0.7 - 0.9 GeV, we obtain 

mi 3(1"^) =0.032 ±0.004 GeV^ . (48) 

2'2 

After combining with (0), we found that the correction to the masses due to the mixing 
is negligibly small. 

In conclusion, we have calculated the l/mq correction to the masses of lowest excited 
heavy mesons from sum rules within the framework of the HQET. This study refines 
the leading order analysis f^. From (42)-(47), together with (pO|), we can obtain the 
corrections to the masses of excited D mesons. For the doublet (0+, 1"*"), our final results 
read 

- (mD* ± ^mD'J = mc ± A ± — [(0.46 ± 0.08) GeV^] , (49) 

mD[ - mo* = — [(0.023 ± 0.003) GeV^] , (50) 
when the currents (|1^) and (p^ ) are used and 

J (mo* ± 3mD[) = me ± A ± — [(0.48 ± 0.10) GeV^] , (51) 

rriD^ - mo* = — [(0.014 ± 0.003) GeV^] , (52) 



m^ 
when the currents (|l^) and ([T6|) are used. As for the doublet (1"'",2"'"), the result is 



- (3mD, + 5mD*) = m^ ± A ± — [(0.50 ± 0.10) GeV^l , (53) 

8 ^ nic 

ruD* -mD, = — [(0.040 ± 0.006) GeV^l , (54) 

^ rric 

where the values of A in three cases are given in [§|. Note that we have reasonably ne- 
glected the renormalization coefficient Cmag of the chromomagnetic operator for charmed 

11 



meson in the discussion. We should also note that the given errors present only the 
uncertainties of T and uo^ in the working regions. 

The above results for the spin splittings are not theoretically reliable because we have 
not included the Feynman diagrams which are of the one loop order in the perturbation 
theory and may be of values comparable to that of the gluon condensate term leading 
in our approximation. Nevertheless, they are in approximate agreement with the ex- 



perimental values for the mass splitting between Di and -Dg 113 ^^ ^^ take rric ^1.3 
GeV. 
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Figure Captions 



Fig. 1. Feynman diagrams contributing to the sum rules for K, in 
the coordinate gauge. The gray square corresponds to the insertion of 
kinetic energy operator at 0{1/Mq) in the HQET Lagrangian. The right 
vertex is put to the origin in coordinate space. There are no figures (c) 



(d) (h) when the currents (|13| ) and (14) are used. 



Fig. 2. Feynman diagrams contributing to the sum rules for S, in the 
coordinate gauge. The gray square corresponds to the insertion of chro- 
momagnetic interaction operator at 0{\/Mq) in the HQET Lagrangian. 

Fig. 3. The sum rules for K as a, functions of the Borel parameter T for 
different values of the continuum threshold ooc- (a) (b) for (0^, 1^) dou- 
blet with the currents ( [l^ (14) and (15) (16) used respectively. Prom top 



to bottom the curves correspond to tOc = 2.9,2.7,2.5 GeV respectively, 
(c) for (1^,2+) doublet. From top to bottom the curves correspond to 
oJc = 3.1,2.9,2.7 GeV respectively. 

Fig. 4. The sum rules for S as a functions of the Borel parameter T for 
different values of the continuum threshold uoc- (a) (b) for (0^, 1^) dou- 
blet with the currents ([T^ (14) and ([l5|) (|l6|) used respectively. From top 



to bottom the curves correspond to uOc = 2.9, 2.7, 2.5 GeV respectively. 
(c) for (l+,2+) doublet. From top to bottom the curves correspond to 
Ldc = 3.1, 2.9, 2.7 GeV respectively. 

Fig. 5. The sum rules for mi 3 (1^) as a functions of the Borel parameter 

2 ' 2 

T for different values of the continuum threshold Wc (from top to bottom: 
Wc = 3.0, 2.8, 2.6). 
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